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Improving the low-lying spectrum of the overlap kernel 
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The action of the overlap-Dirac operator on a vector is typically implemented indirectly through a multi-shift 
conjugate gradient solver. The compute-time this takes to evaluate depends upon the condition number k of the 
matrix that is used as the overlap kernel. We examine the low-lying spectra of various candidate kernels in an 
effort to optimise ft, thereby speeding up the overlap evaluation. 


1. Introduction 

Overlap fermions are a realisation of chiral 
symmetry on the lattice. Given some reasonable 
Hermitian-Dirac operator H, we can deform H 
into a chiral action through the overlap formal¬ 
ism, 

Do = \{l + l^e{H)), e{H) = ^ ( 1 ) 

Unfortunately, the matrix sign function e{H) is 
difficult to evaluate and is typically approximated 
by a sum over poles which can be evaluated 
using a multi-shift conjugate gradient (CG) solver 

. This is an iterative approximation where the 
number of iterations for a given accuracy depends 
upon the condition number of the kernel, k{H) = 
max /A min I ■ 

Usually the Hermitian Wilson-Dirac operator 
is used as the overlap kernel. Its low-lying spec¬ 
trum is characterised by a handful of isolated 
eigenmodes which can be very small, increasing 
K unacceptably. These eigenmodes can be pro¬ 
jected out of the basic operator, reducing its con¬ 
dition number to a numerically acceptable level, 
and then dealt with explicitly]^. Unfortunately, 
as the spectrum rapidly becomes dense, project¬ 
ing out low-lying modes can only help one so far. 

An alternative is to use a kernel with an im¬ 
proved spectrum, that is, where the region of 
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dense modes is shifted away from zero. In this 
paper we consider six different actions as poten¬ 
tial kernels. The first three actions are the stan¬ 
dard Wilson with and without the addition of the 
clover term at both tree-level and with a mean- 
field improved coefficient. The last three actions 
are derived from the first three actions by APE- 
smearing the irrelevant terms. The details of the 
actions are given in [Q, in particular, the FLIG 
(Fat Link Irrelevant Clover) ||] action is defined 
there. 

2. Spectral Flow Comparison 

In order to test the merits of each of our 
proposed actions, we calculate the spectral flow 
of each of them to see if our reasoning regard¬ 
ing their low-lying spectra is correct. From the 
quadratic form of the lower bounds as a function 
of TO, and based upon results given in [|j, we ex¬ 
pect there to be some peak value of to for which 
the gap around zero is the largest. We calculated 
the flow of the lowest 15 eigenvalues as a function 
of TO for two ensembles oflOL = 8^xl6 mean- 
field improved Symanzik configurations, one fine 
ensemble at [3 = 4.80 and one coarse ensemble at 
[3 = 4.38. The following flow graphs allow us to 
see the to value for the biggest gap, and also al¬ 
low us to compare the different actions. As we are 
interested in the magnitude of the low-lying val¬ 
ues rather than their sign, we plot |A| vs to. The 
small physical extent of the fine lattices leads to 
qualitatively different behaviour in the spectral 
flow from the coarse lattice. 
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Figure 1. Spectral flow of the Wilson action at 
(3 = 4.80 (top) and 4.38 (bottom). 


Figure 2. Spectral flow of the clover action at 
f3 = 4.80 (top) and 4.38 (bottom). 


We begin by examining the flow of the Wilson 
action, in Figure |^. We see on the fine lattices 
that there is a reasonable gap away from zero and 
a handful of isolated low-lying eigenmodes before 
the spectra becomes dense. On the coarse lat¬ 
tices the Wilson spectra is very poor, with a high 
density of very small eigenmodes. 

The addition of the clover term (at Csw = 1) 
provides some improvement (Figure |^), shifting 
the flow upwards on the fine lattices and moving 
the peak gap to m = 1 as expected. The clover 
action still performs poorly on the coarse lattices 
however, as the presence of many small eigen¬ 
modes persists, although their density is clearly 


reduced. 

Mean held improvement assists the basic clover 
action somewhat (Figure spreading the spec¬ 
trum upwards, although the lowest modes are 
not raised significantly. The mass value at which 
the peak gap occurs has moved significantly from 
TO = 1 to m = 0.6. As mentioned earlier, es¬ 
sentially all MFI does in this case is to change 
the value of Csw to 1.0/mq, pushing it towards its 
non-perturbative value. 

Modifying the Wilson action by smearing the 
irrelevant operators provides a considerable im¬ 
provement (Figure]^), enhancing the gap around 
zero on the fine lattices significantly. While there 




























3 



Figure 3. Spectral flow of the MFI clover action Figure 4. Spectral flow of the fat Wilson action 
at /3 = 4.80 (top) and 4.38 (bottom). at /3 = 4.80 (top) and 4.38 (bottom). 


are still some small modes present on the coarse 
lattices, their density has been greatly reduced, 
and the spectral flow now has a clear division be¬ 
tween the isolated low-lying modes and the modes 
where the spectral density becomes high. The lat¬ 
ter are well separated from zero. Smearing was 
performed with a = 0.7 and riape = 12 smearing 
sweeps. 

The spectral flow of the fat clover action clearly 
demonstrates the superiority of clover-improved 
actions (Figure The gap around zero is en¬ 
hanced again over the fat Wilson action, and the 
number of isolated low-lying modes is reduced. 
The effect is very pronounced on the coarse lat¬ 


tice. 

As the fat links are already close to unity, the 
addition of mean field improvement only affects 
the fat clover flow slightly (Figure ||), raising the 
gap around zero a little and spreading the eigen¬ 
values upwards slightly also. The low-lying den¬ 
sity is again good in this case and far superior to 
that of the Wilson in Figure |^. On the fine lat¬ 
tices we note that we have a peak gap at m = 1.4 
of 0.5! On the coarse lattices the region where 
the eigenmodes start becoming dense occcurs at 
about 0.35 at the same mass value, indicating 
that the FLIC action should be extremely well- 
conditioned on both ensembles, especially so on 
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Figure 5. Spectral flow of the fat clover action 
at /3 = 4.80 (top) and 4.38 (bottom). 


the coarse lattice if we project out the isolated 
low-lying modes. 

Additionally, we tested the dependence of the 
fat clover and FLIC actions upon the amount of 
smearing done (Figures |^,^. As stated in 0, we 
only effectively need to vary the product ariape, 
so we fix a at 0.7 and vary riape between 0 and 12. 
We observe that in both the case of the fat clover 
and FLIC actions the initial 4-6 sweeps have a 
significant effect, but past 6 sweeps the effect is 
marginal, with the low lying density remaining 
constant and the eigenvalues being compressed 
very slightly downwards. 


Figure 6. Spectral flow of the FLIC action at 
f3 = 4.80 (top) and 4.38 (bottom). 


3. Conclusion 

The addition of the clover-term at tree-level or 
with a mean field improved coefficient does not 
drastically improve the low-lying spectrum of the 
Wilson action at an optimal choice of m. How¬ 
ever, APE smearing the irrelevant-terms in the 
Wilson action both with and without the clover 
term provides a marked improvement in the low- 
lying spectrum, reducing the number of isolated 
low-lying modes and shifting the spectrum away 
from zero. The combination of APE-smearing 
and the clover term at a mean-field improved co¬ 
efficient provides the best overall low-lying spec- 
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Figure 7. Dependence of the fat clover spec¬ 
trum at /3 = 4.80 (top) and 4.38 (bottom) on the 
number of APE smearing sweeps. 


tra. Thus we conclude that the FLIC action is the 
best candidate for the overlap kernel out of the 
six actions tested. These results are built upon in 
m where a comprehensive numerical comparsion 
is undertaken. 
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